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Abstract. Recently, nonlocality and Bell inequalities have been used to investigate 

i-G . quantum phase transitions (QPTs) in low-dimensional quantum systems. Nonlocality 

can be detected by the Bell-CHSH function (BCF). In this work, we extend the study 

of BCF to the QPTs in matrix product systems (MPSs). In this kind of QPTs, the 

C^ , ground-state energy keeps analytical in the vicinity of the QPT points, and is usually 

^ . ' called the MPS-QPTs. For several typical models, our results show that BCF can 

signal the MPS-QPTs very well. In addition, we find BCF can capture signal of 
QPTs in unentangled states and classical states, for which other measures of quantum 
correlation (quantum entanglement and quantum discord) fail. Furthermore, we find 

O^ ■ that in these MPSs, there exists some kind of quantum correlation which cannot be 

VO ' characterized by entanglement, or by nonlocality. 
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1. Introduction 

Quantum phase transition (QPT) is a very interesting phenomenon in many- 
body quantum systems. [Tj Compared with the classical phase transitions driven by 
thermodynamic fluctuation, QPTs occur at zero temperature, thus the thermodynamic 
fluctuation is absent. In fact, QPTs are driven by the so-called quantum fluctuation. 
For a quantum system described by a Hamiltonian H{g) with g the tuning parameter, 
the ground-state property (i.e. the ground-state energy) of H{g) may show qualitative 
change at some point gc, then a QPT occurs. Most QPTs, such as various magnetization 
transitions in spin models, [HE] can be investigated by traditional order parameters. It 
needs mention that, some exotic QPTs, such as the topological QPT, [3] cannot be 
described by local order parameters. 

In the vicinity of the QPT point, long-range correlations would develop in the 
ground state. Thus it is expected that quantum correlation plays a central role in 
the QPTs. [TJ El El [H m El [2] Quantum entanglement is the most famous measure of 
quantum correlation. For various models, it has been found that the entanglement is 
singular in the vicinity of the QPT points, which is usually related to the singularity 
of the ground-state energy at the QPTs. [2l [8] Nonlocality is another aspect of quantum 
correlation, and can be indicated by the violation of Bell inequalities, such as the 
famous Bell-Clauser-Horne-Shimony-Holt (CHSH) inequality. [IDl El El IIIl [iSl El [I3] 
For a long time, nonlocality and entanglement were regarded as identical. In fact, for 
pure two-qubit states, it has been proved rigorously by Gisin that the two are indeed 
equivalent. [14] However, for mixed states, it's found that an entangled state may not 
violate any Bell inequality, i.e., an entangled state is not necessarily a non-local state. [T5] 
Thus, nonlocality and entanglement turn out to be two different aspects of quantum 
correlation. 

Quite recently, it has been found that the Bell-CHSH function (BCF), which is 
associated with the Bell-CHSH inequality, can serve as a useful QPT detector, even 
for topological QPT[3l[l6] and Kosterlitz-Thouless QPT[5]. It should point out that a 
complete understanding about the features of BCF in detecting QPTs has not yet been 
reached. 

In this work, we make a further step by considering an exotic type of QPTs occurring 
in the so-called matrix product states(MPSs).[23l m [211 |22l [HI |19l [20] In this kind of 
QPTs, the ground-state energy remains analytic in the entire parameter space, thus 
the situation is different from traditional QPTs.[T9] In order to distinguish them from 
traditional QPTs, we usually call them MPS-QPTs. The MPSs have always provided 
a valuable test-bed for understanding the features of various QPT detectors, including 
quantum entanglement . [T71 [22] quantum discord [23] and quantum fldelity[2Tl[22]. In this 
work, we investigate the ability of BCF to detect QPTs by considering MPS models. 
Firstly, as we will show, several typical MPSs display clearly the features of BCF in 
detecting QPTs. Secondly, MPSs help us understand the role of quantum correlation, 
quantum entanglement, and nonlocality in QPTs. 
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2. Bell-CHSH inequality, entanglement concurrence, and quantum discord 

Bell-CHSH inequality is the simplest nontrivial Bell inequality. First, let's define the 
CHSH operator as B = Ai Bi + Ai B2 + A2 iS) Bi - A2 B2, where A = «i ■ <? 
and Bi = bi ■ a, with Si and bi unit vectors and a = {a^, (Ty, cTz). Then for any realistic 
and local two-qubit state p2, the Bell-CHSH inequality reads \{B)\ = |Tr(p2i?)| < 2. 
1(5)1 depends upon the vectors Oj and bi, and one can optimize |(-B)| over all vectors 
Si and bi to get the maximum value B{p2) = maxr^. gi |(i?)|.[12l EOl E] We will refer to 
B[p2) as the Bell-CHSH function (BCF) in this work. For some state p2, if it turns out 
that B{p2) > 2, we usually say that the Bell-CHSH inequality is violated, which means 
that the state p2 cannot be described by a realistic local theory, in other words, it is 
non-local. 

Alternatively, it's found by Horodeckis that for any two-qubit state 
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BCF can be expressed by a closed analytical formula. |13] For convenience, here we use 
two works X and o to denote the elements of p2, and for all the models considered in 
this work, it holds that Oij = 0. To calculate BCF, one first defines a 3 x 3 matrix L as 
Lij{p2) = Tr[/92 ■ (3"j (g) (Tj], with {ai, (T2, ds} the Pauli matrices. Then the BCF is given 
by B{p2) = 2y/u + V, with u and v the two largest eigenvalues of the symmetric matrix 

In order to understand the features of nonlocality, we will compare it to two 
closely related measures for bipartite correlation, i.e., entanglement concurrence |27] and 
quantum discord[2ll [25]. Concurrence describes the entanglement between two spins. 
Let's denote p2 as the spin-flipped matrix for the two-qubit density matrix p2, i.e., 
P2 = (^y® ^ypy^y ® cfyi then the concurrence is given by C = max{0, pi — P2 — fJ-s — /U4}, 
where pi are the square roots of the eigenvalues of P2P2 in decreasing order. For 
separable states (in other words, unentangled states), the concurrence would vanish, 
and for maximum entangled states, the concurrence is 1. 

Nonlocality and entanglement are two aspects of quantum correlation. Recently, 
quantum discord is proposed to characterize all the quantum correlation present in the 
system. [231125] Its definition is based on two quantum versions of the classical correlation. 
For a classical system AB composed of two subsystems A and B, the total correlation 
can be expressed as Ia,b = Ha + Hb — Hab, or alternatively, Ja,b = Ha — Ha\b, with 
Ha, Hb and Hab the Shannon entropy, and Ha\b the conditional entropy. Ia,b and 
Ja,b are equal to each other, however, their quantum versions are found to be non- 
equivalent from each other, and the difference is used to define the quantum discord. 
By replacing the Shannon entropy and the conditional entropy with the von Neumann 
entropy and quantum conditional entropy, respectively, Ia,b becomes the quantum 
mutual information Z{pab) = S{pa) + S{pb) — S{pab), and the quantum extension of 
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Ja,b turns out to be the classical correlation J{pab) = maxr^ -,{S{pa) — S{p\{Bk\)}, 
where {Bk} is just a complete set of projectors. |25] Finally, discord is just defined as the 
difference between X{pab) and J^{pab), i-e., D[pab) = ^{pab) — J{.Pab)- For a state 
containing quantum correlation, discord is generally non-zero, while for classical states, 
X{pab) and J{pab) would just reduce to Ia,b and Ja,Bi respectively, thus the discord 
vanishes. 

From the above descriptions, one can see that the entanglement concurrence 
and quantum discord would simply be zero in separable states and classical states, 
respectively. pBJ As we will show in the next section, BCF can capture the signal of 
QPTs in these two situations. 

3. QPTs in MPSs 

In this section, we firstly give a brief introduction to MPSs, then we investigate the 
BCF and nonlocality at MPS-QPTs in several typical models. 

An MPS containing A^ sites (or cells) is defined in the following matrix product 
form [19] 

|^(<7))= Y. Tr(A,...A^)Ki,...,z;v), (2) 

where Ai are D x D matrices, j = 1, ..., N labels the sites, and ij = 1, ...,d denoting the 
degree of freedom for site j. The matrices A^ := Ai{g) depend on the parameter g. 

Usually, for low-dimensional quantum systems, it is difficult to express the ground- 
state wavefunction exactly in an explicit form. However, for an MPS \ilj{g)), one 
can construct a parent Hamiltonian H{g), which guarantees the state {ipig)) be the 
ground state of H{g)- More importantly, in the thermodynamic limit A^ — )■ oo, as the 
change of g, the ground state lipig)) may undergo a novel type of transition, such that 
local observables are singular and the correlation length is divergent, with the ground- 
state energy keeping analytic (In traditional QPTs, the ground-state energy would be 
singular. )[T9]. We usually say that the system undergoes an MPS-QPT. 

For a given MPS, the reduced density matrix of any subsystem in the system can 
be obtained with the help of transfer matrix technique. |20l [TSl [T9] First, let's define the 
transfer matrix E as 

E = Y^A:0A, (3) 

4 = 1 

then the reduced density matrix of k adjacent sites is given by 

Pji,...,ifcji,...,jfc(A) = ( J7N\ ■ ^ -^ 

One can further prove that, for two-site correlation, the correlation length is given 
by ^ = i^rx^i^ p where Ai and A2 are the first and the second largest eigenvalue of the 
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Figure 1. (Color Online) The BCF (B), concurrence (C), and discord (D) as a function 
of a; = tS-^ for the runs; of the ladder. The concurrence has been investigated in Ref. 
|18) and is shown here just for comparison purpose. 



transfer matrix E. Any level crossing between Ai and A2 indicates a divergent correlation 
length, in other words, an MPS-QPT.|19] 

In this work, we only deal with bipartite correlations. For several typical MPS 
models, we use Eq. (^ to obtain the reduced density matrix for the concerned two-site 
subsystem of the models. Then we determine the BCF, and research the behavior of 
BCF, concurrence and discord in the MPS-QPTs. 



3.1. Spin ladder with four-body interaction 

As the first example, we consider an MPS with Ai^2 
/ ' 





A 



1 



. \TE] Its parent Hamiltonian describes a spin s = | ladder model with 



SO (2) symmetry. Every rung of the ladder contains two spins with d = 4, thus four 
matrices are used to define the MPS. The system contains two-body bond interactions 
and four-body plaquette interactions. [H] As the Hamiltonian is too long, we would not 
show it in this work. It has been found that the ladder undergoes an MPS-QPT at 
g = 0, where the spin-spin correlation function of the ladder shows a singularity and 
the largest two eigenvalues of the transfer matrix have a level crossing. 

Now we try to use quantum correlations to find the signal for the MPS-QPT. Let's 
just consider a single rung in the ladder. The elements of the reduced density matrix 
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P2 turn out to be 

{x22,a;23, 2^32, 3:33} = 1 , (5) 

{0,2:14^X41} = 

where ^ '■= ^- We have determined BCF on the rung and shown it in Fig. [T] The 

BCF shows a singularity at the MPS-QPT point (? = 0, thus it can be used to detect 

the MPS-QPT in this model. 

In addition, the first-order derivative of BCF is discontinuous at a; = ±2. Detailed 

analysis shows that the singular points aX x = ±2 are due to the mathematical definition 

/vfrunff) 
of BCF, rather than the singularity in p2 . Explicitly, the non-physical singularity 

is induced by the max function in the definition of BCF. In order to calculate BCF, one 

has to find the two largest eigenvalues of the symmetric matrix L^ L. In this procedure, 

a mathematical singularity may emerge. In fact, a max/min function is also involved in 

the definition of the discord/concurrence. As a result, the concurrence and discord can 

also show a non-physical singularity, just as shown in Fig. 1. 

Now let's discuss the feature of quantum correlation in the rung. As indicated by 

the discord in Fig. [H quantum correlation exists for any finite x. In the vicinity of 

the MPS-QPT point, that is, for \x\ < 0.41, it's found that the concurrence is non-zero 

and B > 2, thus the quantum correlation is in the form of both entanglement and 

nonlocality. While for 0.41 < \x\ < 1, the concurrence is non-zero and B < 2, thus the 

quantum correlation is in the form of entanglement without nonlocality. For |a;| > 1, it 

is present neither in the form of entanglement nor in the form of nonlocality. It shows 

clearly that quantum correlation can be manifested by various forms. 

3.2. XYZ interaction model 

Let's consider an MPS with Ai = and Ao = { . After the standard 

\ 1 1 / 1-11/ 



procedure, one can construct its parent Hamiltonian as [20] 

N 

H = Y: J.KK'' + Jy^W + JzK^' - Ba:, (6) 

with J^ = — J + |(1 + g"^), Jy = —J + g, Jz = —J — g and B = 1 — g"^. It is just an XYZ 
interaction chain. It has been proved that the system has an MPS-QPT aX g = 0.[20] 

Firstly, we consider the correlation between two nearest-neighboring spins i and 
i + 1. For (? > 0, the corresponding reduced density matrix p2 is given by 
{xii,a;44} = 5(2 + 651 + 1 

{X22, a;23, 2^32, 3^33, 2^14, 2:41} = {g - if , (7) 

0=1- g^ 

while for g > 0, one finds that 

{a;i4,a;4i} = 5^ + 65 + 1 

{x22, a;23, a;32, a^33, Xu,Xu} = {g- 1)^ ■ (8) 

o = l-g' 
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Figure 2. (Color Online) The BCF (B) as a function of g for any two spins in the 
XYZ interaction model. The concurrence (C) has been investigated in Ref. [3D] and 
is shown here just for comparison purpose. 



Previous studies show that pg* is separable for any g, thus the bipartite 
entanglement between i and z + 1 vanishes and cannot detect the MPS-QPT of the 
system. [20| 123] From Fig. |2] we find that BCF is generally non-zero in the whole 
parameter space and shows a singularity at the MPS-QPT point g = Q. 

In fact, for any two-site subsystems of the model, the reduced density matrix pg '* 
does not depend on the distance r at all. [20] As a result, the concurrence for any two-site 
subsystem is zero, thus bipartite entanglement cannot detect the QPT of the system 
while the BCF can do the job very well. This example shows clearly that BCF can 
detect QPTs in separable states while entanglement fails. 

We observe that the Bell-CHSH inequality is not violated in the MPS-QPT in 
this XYZ interaction model. We'd like to mention that the quantum correlation indeed 
exists in the QPT region, indicated by the discord. [23] It is interesting that the quantum 



correlation present in pg 
nonlocality. 



(i,i+r) 



is neither in the form of entanglement nor in the form of 



3.3. Three-body interaction model 



We consider an MPS with Ai 



and ^2 = 

^ 1 1 y ^ V 

describes a three-body interaction model as follows [19] 

N 



H = Y.J^^z(^x^^ 



-|- J zO ,U ^ 



i=l 



Its parent Hamiltonian 



(9) 



with J^ = {g - 1)2, J^ 
MPS-QPT at £/ = 0.[19j 



2{g'^ — 1), and B = (1 -|- ^f)^. The system undergoes an 
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Figure 3. (Color Online) The BCF (B) as a function of g for nearest-neighboring two 
sites in the three-body interaction model. The concurrence (C) has been investigated 
in Ref. 23 and is shown just for comparison purpose. 



We consider the two-qubit states p2 ^ in the chain with different r. When r 
the reduced density matrix for g > is given by 



{Xii,X44} = ^, {x22,X3-i} ^ 3-^ 



{X23,X32} — -^, {Xi4,X4i} — -^ 
= 9 



(10) 



while ioi g < 0, P2' is a diagonal matrix, with the diagonal entries given by 

{Xii,2;44} = 1, {X22,X33} = -g , (11) 

which denotes a classical state. The BCF is shown in Fig. |3l BCF shows a singularity 
at g = 0, thus it can be used to detect the MPS-QPT of the model. In addition, BCF 
is singular at ^f = — 1. Detailed analysis shows that this singular point results from the 
mathematical definition of BCF, rather than a transition in \ip{g)). For g > 0,we observe 
that P2 * never violate the Bell inequality, despite being entangled. In other words, 
the quantum correlation in nearest-neighboring sites is in the form of entanglement, 
rather than nonlocality. 

Next, we consider two sites i and i + r with r > 2. For 5^ > 0, the elements of the 
reduced density matrix /)2 '* ^ are given by 



1 + 



{xu, 2:44} 

{X22, 2:33} 

{Xi4, a;23, 3^32, 3^41/ 

o 



1- fi^l 



4.9 

(1+S)2 
{i,i+r) 



16.9^ 



(12) 



while for (7 < 0, P2 
by 

{Xu, X44} 



is reduced to a diagonal matrix, with the diagonal elements given 

{X22,XSS} = I- Ct^J . (13) 
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Figure 4. (Color Online) The quantum discord (D) and the BCF (B) as a function 
of g for two sites i and i -f- r in the three-body interaction model. 



We numerically found that the concurrence vanishes in the whole parameter space 
for r > 2, which means that p2 '^ is separable. Previous study shows that quantum 
discord may be able to capture the signal of QPT in separable states[26l |23]. We have 
calculated the discord for different r and the result is shown in Fig. Hl^a). One can see 
that the value of discord is very small even for r = 2. As the increase of the distance 
r, the discord decreases rapidly. Finally, for a large r, the two-qubit state /92 '* ^ would 
become a classical state without any quantum correlation, thus neither the concurrence 
nor the discord can signal the MPS-QPT of the system. Then let's study the BCF for 
^(1,1 r) Pj,qj^ pjg gl^b) one sees clearly that BCF shows a singularity at the QPT point 
g = for any finite r. Thus, BCF is able to capture the singularity in classical states, 
for which both the concurrence and the discord fail. In addition, our results show that 
the Bell-CHSH inequality is never violated, thus the quantum correlation between the 
non- nearest-neighbor spins, if exists, is neither in the form of entanglement, nor in the 
form of nonlocality. 



4. Summaries and discussions 

In this work, for several typical models, we find that BCF can be used to detect the 
MPS-QPTs very well. The underlying mechanism is as follows. Our discussion applies 
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to MPS-QPTs and traditional QPTs. When a QPT occur, some local observables, 
such as spin-spin correlation functions, would be singular at the QPT point. [28] These 
correlation functions can be used to construct the reduced density matrix p2(fl')-[29] As 
a result, the QPT and the singularity in P2{g) are closely related to each other. The 
BCF, concurrence and discord, defined based upon p2{g), may capture the singularity 
in P2{g), thus, all the three quantities can be used to detect the QPTs. 

In Ref. |5] a general argument for why BCF should be as good as entanglement 
to signal QPTs has already been made. Now let's discuss the advantage of BCF in 
detecting QPTs. BCF can capture singularity in unentangled states and classical states, 
for which entanglement and discord fails, respectively. Quantum discord is defined on 
the classical- quantum paradigm from a measurement perspective, thus it captures all 
the quantumness of correlation in the system. For classical states, discord is zero thus 
loses it function as a QPT indicator. This is just the situation for pg '* ^ with r > 1 
in the three-body interaction model. On the other hand, quantum entanglement is 
defined on the separability-entanglement paradigm. For quantum separable states, it 
is zero thus would lose the signal of QPT in these states. This is just what happens 
in any two-site subsystem in the XYZ model. BCF is used to detect the nonlocality 
of a state. For non-local states, it would be larger than 2, while for most local states, 
it does not vanish. In extreme cases, even if the state is a classical state thus p2 is 
diagonal, the BCF would turn out to be 5 = 2|a:;ii -|- 0:44 — X22 — a^ssl, which is still 
non-zero unless Xu + 2:44 = X22 + a^ss- As a result, BCF can detect the singularity in 
various density matrices, including separable/entangled states, local/non-local states, 
and classical/quantum states. 

Now we clarify the drawback of BCF. First of all, for a general state, the calculation 
of BCF is difficult, which greatly limits its application. In addition, as shown in x = ±2 
of Fig. [1] and (7 = — 1 of Fig. [3|, the mathematical definition of BCF can introduce 
non-physical singularity, which is not related to the singularity in p{g). As a result, the 
singularity of BCF just can be used to detect, rather than to determine, an MPS-QPT. 
However, we'd like to mention that both the concurrence and the discord have similar 
disadvantage (see Fig. [T]as examples), which has already been discussed in some other 
studies. |26l[H] 

Let's discuss the form of quantum correlation in the MPS-QPTs of these models. 
In this work, we have only dealt with bipartite correlations. We use the discord to 
identify the existence of quantum correlation, and then describe the nature of quantum 
correlation through the analysis of entanglement and nonlocality. In the vicinity of 
MPS-QPT point in the ladder model, the quantum correlation in p2 is in the 

form of both entanglement and nonlocality. In the three-body interaction model, for 
P2'* , the quantum correlation is in the form of entanglement without nonlocality, 
while for p2 * ^ with r > 1, the quantum correlation, if exist, is neither in the form of 
entanglement nor nonlocality. In the XYZ interaction model, for any two spins in the 
chain, the quantum correlation is neither in the form of entanglement nor nonlocality. 
From one hand, our results show that when MPS-QPTs occur, the two-site quantum 
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correlation can show very rich nature combined with entanglement and nonlocality. 
From another hand, it reveals that entanglement and nonlocality are not the only aspects 
of quantum correlation, and there exists some kind of quantum correlation which cannot 
be characterized by entanglement, nor by nonlocality. 

Finally, as we have shown in this work, the Bell-CHSH inequality is violated just in 
the QPT of the four-body interaction ladder model. In fact, as far as we know, in all the 
previous works, p^ [5l |6] when QPTs occur in infinite models, the density matrices of 
two-qubit subsystems never violate the Bell-CHSH inequality. Thus, the ladder model 
reported in this work may be the first one to present such a behavior. In QPTs in 
many-body systems, as an important aspect of quantum correlation, nonlocality should 
play a central role, however, it turns out that bipartite nonlocality is not a common 
form of quantum correlation present in these one- dimensional systems. 

We'd like to mention that, for many-body systems, which are naturally multipartite, 
it would be more natural for quantum correlation to present in the form of multipartite 
nonlocality. [351 l30l l3T] It has been found that Bell inequalities can be used to test 
multipartite nonlocality. [32] In addition, effective approaches to detect or even quantify 
multipartite nonlocality have been proposed. [331 El] 

However, the relevance of multipartite nonlocality in QPTs remains unknown. As 
bipartite nonlocality is not favored at the QPT points in many one- dimensional models, 
it would be interesting to clarify whether multipartite nonlocality is significant in QPTs. 
Considering the simple product form of MPSs, we believe MPS models would be very 
useful to study this important issue. 

After finishing this manuscript, we become aware of a related work [36j by Oliveira 
et. al. The authors have proposed a general explanation for why the Bell inequality 
is not violated in most translation invariant systems. Then they have shown that the 
inequality can be violated for models with translation symmetry breaking. Our results 
of the spin ladder model suggest that, for two spins located in a unit-cell of a complex 
lattice, the Bell inequality can still be violated in a translation invariant system. 
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